BUNDESREPUBLIK DEUTSCHLAND 




Prioritatsbescheinigung uber die Einreichung 
einer Patentanmeldung 



,ffv Aktenzeichen: 102 50 922.0 

Anmeldetag: 31 . Oktober 2002 



Anmelder/lnhaber: Siemens Aktiengesellschaft, Munchen/DE 



Bezeichnung: Verfahren zur Ermittlung des ADC-Koeffizienten in 

der Diffusions-gewichteten Magnetresonanz-Bild- 
gebung bei Verwendung von Steady-State-Sequen- 
zen 



IPC: G 01 R 33/563 

Die angehefteten Stucke sind eine richtige und genaue Wiedergabe der ur- 
sprunglichen Unterlagen dieser Patentanmeldung. 




A 9161 

06/00 
EDV-L 



200204698 



Beschreibung 



Verfahren zur Ermittlung des ADOKoef f izienten in der Diffu- 
sions-gewichteten Magnetresonanz-Bildgebung bei Verwendung 
5 von Steady-State-Sequenzen 

Die vorliegende Erfindung bezieht sich allgemein auf die 
Kernspintomographie (Synonym: Magnetresonanztomographie, MRT) 
wie sie in der Medizin zur Untersuchung von Patienten Anwen- 
10 dung findet. Dabei bezieht sich die vorliegende Erfindung 
insbesondere auf ein Verfahren zur Ermittlung des ADC- 
Koeff izienten in der Dif fusions-gewichteten Magnetresonanz- 
^ Bildgebung. 

15 Die MRT basiert auf dem physikalischen Phanomen der Kernspin- 
resonanz und wird als bildgebendes Verfahren seit uber 15 
Jahren in der Medizin und in der Biophysik erfolgreich einge- 
setzt. Bei dieser Untersuchungsmethode wird das Objekt einem 
starken, konstantem Magnetfeld ausgesetzt. Dadurch richten 

20 sich die Kernspins der Atome in dem Objekt, welche vorher re- 
gellos orientiert waren, aus . Hochf requenzwellen konnen nun 
diese „geordneten* Kernspins zu einer bestimmten Schwingung 
anregen. Diese Schwingung erzeugt in der MRT das eigentliche 
Messsignal, welches mittels geeigneter Empf angsspulen aufge- 

25 nommen wird. Durch den Einsatz inhomogener Magnetf elder , er- 
^ zeugt durch Gradientenspulen, kann dabei das Messobjekt in 

alle drei Raumrichtungen raumlich kodiert werden was im All- 
gemeinen als „Ortskodierung* bezeichnet wird. 

30 Bei der Beurteilung pathophysiologischer Vorgange, insbeson- 
dere im menschlichen Gehirn, beispielsweise bei einem Schlag- 
anfall, erweist sich eine verhaltnismaftig neue MR-Technik als 
besonders wirksam: die Dif f usions-gewichtete Magnetresonanz- 
Tomographie . 

35 

Diffusion entsteht durch die thermische Translationsbewegung 
von Molekulen. Dabei handelt es sich urn einen Zuf allsprozefi , 
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der auch als Brownsche Molekularbewegung bezeichnet wird. Die 
zuruckgelegten Distanzen der in dif f usionsgewichteten MRT- 
Messungen betrachteten Molekule sind sehr gering; z.B. dif- 
fundieren Wassermolekule uneingeschrankt typischerweise in 
5 jede beliebige Richtung tiber eine Distanz von ca. 20|nm in 

100ms bzw. 60[xm in Is. Diese Distanzen liegen in der Grofien- 
ordnung einzelner Zellen insbesondere des menschlichen Zell- 
gewebes . Durch den Einsatz ausgesprochen starker Magnet- 
Gradientenf elder (sogenannter Dif fusionsgradienten) , die bei 

10 dieser Technik zusatzlich zu den oben genannten ortskodieren- 
den Gradientenf eldern permanent oder auch pulsweise angelegt 
werden, macht sich eine kollektive Dif f usionsbewegung der je- 
t 'X weiligen Molekule (insbesondere Wasser) in einer Abschwachung 
des Magnetresonanzsignals bemerkbar. Bereiche in denen Diffu- 

15 sion auftritt werden daher mehr oder weniger als dunkle Be- 
reiche im tatsachlichen MRT-Bild gekennzeichnet , je nach 
Starke der Diffusion, welche von verschiedenen Faktoren ab- 
hangt. Die genaue Theorie der Signalentstehung in der Diffu- 
sions-gewichteten Bildgebung wird in der Figurenbeschreibung 

20 anhand von Zeichnungen erlautert. 

Ein Problem in den Anfangen der Dif f usions-gewichteten Bild- 
gebung war die ausgesprochene Empf indlichkeit auf nicht- 
dif fusionsartige Bewegungsarten wie: Herzbewegung, Atembewe- 

25 gung, usw., und die damit verbundenen Bewegungen wie bei- 

spielsweise der Hirnpulsation (Bewegung des Hirns im Liquor) . 
Den Einsatz der Dif f usionsbildgebung als klinische Untersu- 
chungsmethode hat erst die stetige Weiterentwicklung schnel- 
ler Mefitechniken, wie beispielsweise der Echoplanaren Bildge- 

30 bung (engl: Echo-Planar-Imaging, EPI) moglich gemacht . EPI 
ist eine ausgesprochen schnelle MeUmethode in der MRT . Bei 
Verwendung von Single-Shot-Echo-Planar-Imaging (SSEPI)- 
Sequenzen lassen sich Bildartef akte, die aufgrund von unver- 
meidbaren Bewegungs typen entstehen, reduzieren bzw. vermei- 

35 den. Bewegungen wie sie bei konventionellen dif f usionsgewich- 
teten Bildgebungssequenzen auftreten konnen mit SSEPI gewis- 
sermaBen "eingef roren" werden. Ein Nachteil, bedingt durch 
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die Art der Phasenkodierung einer SSEPI-Sequenz ist jedoch 
die sehr starke T 2 *- (T 2 * ist die Zerf allsdauer der Quermagne- 
tisierung unter Berucksichtigung von lokalen Magnet feldinho- 
mogenitaten) bzw die sehr starke Phasenempf indlichkeit . Bei- 
5 des resultiert in starken Bildausloschungs- bzw. Verzerrungs- 
artefakten insbesondere bei der Korperbildgebung mit typisch 
kurzen T 2 -Zeiten des menschlichen Gewebes. 



"Nicht-EPI-Sequenzen" (allgemein als Steady-State-Sequenzen 
10 bezeichnet), wie beispielsweise FISP (engl.: Fast-Imaging- 
with-Steady-Precession) und PSIF (Umkehrung von FISP) , ver- 
wenden generell die tibliche spin warp-Phasenkodiertechnik 
(diskrete Phaseninkrementierung durch einen Phasenkodiergra- 
dient) und sind bzgl. der oben beschriebenen Artefakte unemp- 
15 findlich. Ublicherweise wird bei einer derartigen Sequenz ein 
monopolarer (positiver oder negativer) in der Regel pulsarti- 
ger Dif f usionsgradient geschaltet mit einem HF-Anregungspuls 
von a<90°. Unabhangig vom Dif fusionsgradienten hat ein der- 
artiger HF-Puls a drei unterschiedliche Einwirkungs- 
20 Eigenschaf ten auf einen Magnetisierungsvektor : 

1. Kippen des Magnetisierungsvektors relativ zur longitudina- 
len Richtung (z-Achse) um den Flipwinkel a, 

2. Invertieren des Magnetisierungsvektors um 180° , und 
25 3. Keinerlei Wirkung auf den Magnetisierungsvektor. 

Wie spater anhand der Figuren ausfuhrlich erlautert wird, er- 
geben sich aufgrund dieser drei Eigenschaf ten des HF- 
Anregungspulses unterschiedlich verzweigte Phasenver lauf e der 

30 Langs- und Quermagnetisierung (auch "Echo-Pf ade" genannt) de- 
ren jeder eine unterschiedliche Dif f usionszeit A ± aufweist. 
Abhangig vom Phasengang bzw. der Phasenhistorie, welche durch 
einen anliegenden Dif fusionsgradienten nochmals eine Auffa- 
cherung erfahrt, ist auch die entsprechende T x - und T 2 -Wich- 

35 tung wahrend des gesamten Phasenverlauf es eines bestimmten 
Echopfades. Insgesamt ergibt sich ein MRT-Signal, welches 
durch das Zusammenlauf en einer Vielzahl von Echopfaden gebil- 
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det wird und daher eine Komplexitat aufweist, die nicht mehr 
berechnet werden kann. 

Ein dif fusionsgewichtetes MRT-Bild ergibt sich aus dem die 
5 Diffusion charakterisierenden pixelweise bestimmten D ADC -Wert 
(engl. : Apparent-Diffusion-Coefficient, ADC-Koef f izient ) der 
sich aus den Melisignalen der jeweiligen verwendeten Sequenz 
sowie aus dem das Experiment charakterisierenden b-Wert be- 
rechnen lafit. Hierzu miissen bei einer dif fusionsgewichteten 

10 PSIF-Sequenz mit unipolarem Dif f usions-Gradientenpuls bei der 
Integration aller Echopfade die jeweiligen Ti- und T 2 -Werte 
sowie die Flipwinkelverteilung <x(z) des verwendeten realen 
^/ HF-Puls uber die selektierte Schicht exakt bekannt sein 

[R.Buxton, J. of Magnetic Resonance in Medicine 29, 235-243 

15 (1993) ] . Dies ist bei den bisher verwendeten (Steady-State-) 
Sequenzen nicht der Fall weshalb eine genaue Berechnung des 
D ADC -Wertes unter Berucksichtigung einer exakten T 2 -Gewichtung 
nicht moglich ist. Die Folge sind dif f usionsgewichtete Bilder 
die starke Artefakte aufweisen. 

20 

Aus diesem Grunde schlagen Y.Zur, E.Bosak, N.Kaplan im Jour- 
nal of Magnetic Resonance in Medicine 37, 716-722 (1997) vor 
anstelle eines monopolaren Dif f usionsgradienten einen bipola- 
ren Dif f usionsgradienten einzusetzen. Hierdurch wird das 
25 durch den Di ff usionsgradienten bewirkte Auffachern der Pha- 
Q senverlaufe kompensiert. Die Dif fusionszeit A ist dadurch 

wohldef iniert und eine ADC-Berechnung aus zwei Dif fusionsmes- 
sungen mit hinsichtlich der Amplitude unterschiedlichen Dif- 
fusionsgradienten naherungsweise moglich. 

30 

Allerdings war bereits bekannt, dass bei Verwendung bipolarer 
Dif fusionsgradienten die berechneten Dif f usionskoef f izienten 
immer noch stark vom jeweiligen T 2 -Wert abhangen, was eben- 
falls zu erheblichen Artefakten fuhrt [S.Ding, H.Trillaud et 
35 al. J. Of Magnetic Resonance in Medicine 34, 586-595 (1995)]. 
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M.H.Cho und C.H.Cho haben bereits in dem Konf erenzband Socie- 
ty of Magnetic Resonance in Medicine, p. 911 Amsterdam (1989) 
theoretisch gezeigt, dass diese T 2 -Wichtung bei Verwendung 
eines bipolaren Dif f usionskoef f izienten eliminiert werden 
5 kann, wenn zu dem gemessenen PSIF-Signal (S~) ein entspre- 
chendes FISP-Signal (S + ) bekannt ist und eine entsprechende 
'Formel angegeben . 

Aufgabe der vorliegenden Erfindung ist es, die far diese For- 
10 mel benotigten Messgroiien besonders vorteilhaft zu messen urn 
eine artef akt f reie Berechnung des D ADC -Wertes zu ermoglichen. 



Diese Aufgabe wird gemail der Erfindung durch die Merkmale der 
unabhangigen Anspruche gelost. Die abhangigen Anspruche bil- 
15 den den zentralen Gedanken der Erfindung in besonders vor- 
teilhafter Weise weiter. 



Erf indungsgemafi wird ein Verfahren zur Ermittlung eines dif- 
fusionsgewichteten Bildes in der dif f usionsgewichteten MRT- 
20 Bildgebung vorgeschlagen welches die folgenden Schritte auf- 
weist: 

a) Messen und Speichern eines nichtdif fusionsgewichteten und 
eines dif fusionsgewichteten Datensatzes mittels einer 

25 DESS-Sequenz (Double Echo Steady State Sequenz), wobei fur 

den nichtdif fusionsgewichteten Datensatz zwei Auslesegra- 

"> dienten nacheinander geschaltet werden und fur den diffu- 

sionsgewichteten Datensatz zwischen zwei Auslesegradienten 
ein bipolarer Dif f usionsgradientenpuls zug geschaltet wird, 

30 

b) Berechnen eines dif fusionsgewichteten MRT-Bildes auf Basis 
des nichtdif fusionsgewichteten und des dif fusionsgewichte- 
ten Datensatzes sowie auf Basis eines die dif f usionsge- 
wichtete Messung charakterisierenden Wertes. 

35 



Der bipolare Dif f usionsgradientenpulszug besteht vorteilhaft- 
erweise aus einem positiven und einem negativen Diffusions- 
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gradientenpuls mit der jeweiligen Amplitude G 0 bzw. -G 0 sowie 
der gleichen Pulsbreite 5. Beide Dif f usionsgradientenpulse 
folgen erf indungsgemaii unmittelbar hintereinander . Diese Art 
der Gradientenschaltung kompensiert auftretende Dephasierung- 
5 seffekte. 



Die Berechnung des dif f usionsgewichteten MRT-Bildes erfolgt 
schliefilich durch Bilden des Quotienten einer Kombination des 
dif fusionsgewichteten und des nichtdiff usionsgewichteten Da- 
10 tensatzes, anschlieliender Logarithmierung des Quotienten so- 
wie Gewichtung mit dem der dif fusionsgewichteten Messung ent- 
sprechenden b-Wertes. 



Erfindungsgemafi erfolgt die Bildgebung des dif f usionsgewich- 
15 teten MRT-Bildes durch pixelweise Darstellung eines ADC- 

Koef f izienten Z)^ c der aus den akquirierten Datensatzen gemaB 
der Forme 1 



20 



30 



ebenfalls pixelweise ermittelt wird, wobei 

S* und Sq den Datensatz der nichtdiff usionsgewichteten Mes- 
sung in Form von FISP-Echo-Signalen und von PSIF-Echo- 
Signalen darstellen sowie 
|P5 So W und S~ ijr den Datensatz der dif fusionsgewichteten Messung 

in Form von FISP-Echo-Signalen und von PSIF-Echo-Signalen 
darstellen, und wobei 

b bip den die dif f usionsgewichtete Messung charakterisierenden 
Wert darstellt der nach der Formel 



35 



durch die Beschaf f enheit des bipolaren Di f f usionsgradienten 
und dem gyromagnetischen Verhaltnis y der betrachteten kern- 
resonanten Spinspezies gegeben ist. 
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Die Gesamtmesszeit kann vorteilhaft dadurch wesentlich redu- 
ziert werden wenn das S*-Signal im Vergleich zum S"-Signal 
erf indungsgemafi mit hoherer Bandbreite gemessen wird. 

Urn eventuell auftretende Artefakte zu vermeiden ist es vor- 
teilhaft, beide Signale S + und S~ mit gleicher groBer Band- 
breite zu messen, wobei das S~-Signal durch einen Multigra- 
dientenechozug mehrmals gemessen wird und liber die gemessenen 
S~-Signale gemittelt wird. 

Vorteilhaf terweise bietet es sich an die Mittelung mit dem 
Quadratsummenverf ahren durchzuf uhren . 

Die Akquirierung der Datensatze S~ iJf , S+ iJf , S~ , £ 0 + kann erf in- 

dungsgemaii ebenso durch ein Pro j ektions-Rekonstruktions- 
Verf ahren erf olgen . 

Ferner wird erf indungsgemail ein Kernspintomographiegerat zur 
Durchfuhrung der erf indungsgemafien Ver f ahrenschritte bean- 
sprucht . 

Weitere Vorteile, Merkmale und Eigenschaf ten der vorliegenden 
Erfindung werden nun anhand von Ausf uhrungsbeispielen bezug- 
nehmend auf die begleitenden Zeichnungen naher erlautert. 

Figur 1 zeigt schematisch ein Kernspintomographiegerat, 

Figur 2a zeigt schematisch den zeitlichen Verlauf der Gra- 

dientenpulsstromfunktionen einer Spin-Echo-Sequenz , 

Figur 2b zeigt schematisch die zeitliche Abtastung der k- 
Matrix bei einer Spin-Echo-Sequenz, 



Figur 3 zeigt schematisch den Phasenverlauf (J) der Quermag- 
netisierung sowie die Signalentstehung in einem 
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Spin-Echo-Experiment unter Einflufi eines konstant 
anliegenden Diffusions gradient en, 

Figur 4 zeigt schematisch den Phasenverlauf <\> der sich aus 
5 Magnetisierungsanteilen zusammenset zenden Quermag- 

netisierung bei Rephasierung und dif f usionsbeding- 
ter Dephasierung in einem Spin-Echo-Experiment, . 

Figur 5 zeigt schematisch eine Stej skal-Tanner-dif f usions- 
10 gewichtete Spin-Echo-Sequenz, 

Figur 6 zeigt schematisch die Struktur einer DESS-Sequenz 
^ sowie die jeweiligen Phasenentwicklung der Magneti- 

sierungsvektoren bei FISP und PSIF. 

15 

Figur 7 zeigt eine erste Ausf uhrungsf orm der er f indungsge- 
maften MRT-Mess-Sequenz mit bipolarem Dif f usionsgra- 
dient und FISP- und PSIF-Messung, 

20 Figur 8 zeigt eine zweite Ausf uhrungsf orm der erf indungsge- 
mafien MRT-Mess-Sequenz mit bipolarem Dif f usionsgra- 
dient, FISP-Messung sowie n-facher Multigradienten- 
echo-Messung einer PSIF-Messung, 

25 Figur 9 zeigt eine dritte Ausf uhrungsf orm der er f indungsge- 
maBen MRT-Mess-Sequenz mit bipolarem Dif f usionsgra- 
\ dient und FISP- und PSIF-Messung unter Verwendung 

des Projektions-Rekonstruktions-Verfahrens . 



30 Figur 1 zeigt eine schematische Darstellung eines Kernspin- 
tomographiegerates zur Erzeugung von Gradientenpulsen gemaft 
der vorliegenden Erfindung. Der Aufbau des Kernspintomogra- 
phiegerates entspricht dabei dem Aufbau eines herkommlichen 
Tomographiegerates . Ein Grundf eldmagnet 1 erzeugt ein zeit- 

35 lich konstantes starkes Magnetfeld zur Polarisation bzw. Aus- 
richtung der Kernspins im Untersuchungsbereich eines Objek- 
tes, wie z.B. eines zu untersuchenden Teils eines menschli- 
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chen Korpers. Die fur die Kernspinresonanzmessung erforderli- 
che hohe Homogenitat des Grundmagnetf eldes ist in einem ku- 
gelformigen Meflvolumen M definiert, in das die zu untersu- 
chenden Teile des menschlichen Korpers eingebracht werden. 
5 Zur Unterstiitzung der Homogeni tatsanf orderungen und insbeson- 
dere zur Eliminierung zeitlich invariabler Einfliisse werden 
an geeigneter Stelle sogenannte Shim-Bleche aus f erromagneti- 
schem Material angebracht . Zeitlich variable Einflusse werden 
durch Shim-Spulen 2 eliminiert, die durch eine Shim- 
10 Stromversorgung 15 angesteuert werden. 

In den Grundf eldmagneten 1 ist ein zylinder f ormiges Gradien- 
^ tenspulensystem 3 eingesetzt, das aus drei Teilwicklungen be- 
steht. Jede Teilwicklung wird von einem Verstarker 14 mit 
15 Strom zur Erzeugung eines linearen Gradientenf eldes in die 
jeweilige Richtung des kartesischen Koordinatensystems ver- 
sorgt. Die erste Teilwicklung des Gradientenf eldsystems 3 er- 
zeugt dabei einen Gradienten G x in x-Richtung, die zweite 
Teilwicklung einen Gradienten G y in y-Richtung und die dritte 
20 Teilwicklung einen Gradienten G z in z-Richtung. Jeder Ver- 
starker 14 umfaBt einen Digi tal-Analog-Wandler , der von einer 
Sequenzsteuerung 18 zum zeitrichtigen Erzeugen von Gradien- 
tenpulsen angesteuert wird. 



m 



25 Innerhalb des Gradientenf eldsystems 3 befindet sich eine 
Hochf requenzantenne 4, die die von einem Hochf requenz- 
leistungsverstarker 30 abgegebenen Hochf requenzpulse in ein 
magnetisches Wechselfeld zur Anregung der Kerne und Ausrich- 
tung der Kernspins des zu untersuchenden Objektes bzw. des zu 
30 untersuchenden Bereiches des Objektes umsetzt. Von der Hoch- 
f requenzantenne 4 wird auch das von den prazedierenden Kerns- 
pins ausgehende Wechselfeld, d.h. in der Regel die von einer 
Pulssequenz aus einem oder mehreren Hochf requenzpulsen und 
einem oder mehreren Gradientenpulsen hervorgeruf enen Kerns- 
35 pinechosignale, in eine Spannung umgesetzt, die uber einen 
Verstarker 7 einem Hochf requenz-Empf angskanal 8 eines Hoch- 
f requenzsystems 22 zugefiihrt wird. Das Hochf requenzsystem 22 
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umfaUt weiterhin einen Sendekanal 9, in dem die Hochf requenz- 
pulse fur die Anregung der magnetischen Kernresonanz erzeugt 
werden. Dabei werden die jeweiligen Hochf requenzpulse auf- 
grund einer vom Anlagenrechner 20 vorgegebenen Pulssequenz in 
5 der Sequenzsteuerung 18 digital als Folge komplexer Zahlen 
dargestellt. Diese Zahlenfolge wird als Real- und als Imagi- 
naranteil uber jeweils einen Eingang 12 einem Digital-Analog- 
Wandler im Hochf requenzsystem 22 und von diesem einem Sende- 
kanal 9 zugefuhrt. Im Sendekanal 9 werden die Pulssequenzen 
10 einem Hochf requenz-Tragersignal aufmoduliert , dessen Basis- 
frequenz der Resonanzf requenz der Kernspins im MeJJvolumen 
entspricht • 



Die Umschaltung von Sende- auf Empf angsbetrieb erfolgt liber 
15 eine Sende-Empf angsweiche 6. Die Hochf requenzantenne 4 

strahlt die Hochf requenzpulse zur Anregung der Kernspins in 
das MeBvolumen M ein und tastet resultierende Echosignale ab . 
Die entsprechend gewonnenen Kernresonanzsignale werden im 
Empf angskanal 8 des Hochf requenzsystems 22 phasenempf indlich 
20 demoduliert und liber einen jeweiligen Analog-Digital-Wandler 
in Realteil und Imaginarteil des MefJsignals umgesetzt. Durch 
einen Bildrechner 17 wird aus den dergestalt gewonnenen Mefi- 
daten ein Bild rekonstruiert . Die Verwaltung der MeJidaten, 
der Bilddaten und der Steuerprogramme erfolgt uber den Anla- 
25 genrechner 20, Auf grund einer Vorgabe mit Steuerprogrammen 
kontrolliert die Sequenzsteuerung 18 die Erzeugung der je- 
weils gewunschten Pulssequenzen und das entsprechende Abtas- 
ten des k-Raumes. Insbesondere steuert die Sequenzsteuerung 
18 dabei das zeitrichtige Schalten der Gradienten, das Aus- 
senden der Hochf requenzpulse mit definierter Phase und Ampli- 
tude sowie den Empfang der Kernresonanzsignale. Die Zeitbasis 
fur das Hochf requenzsystem 22 und die Sequenzsteuerung 18 
wird von einem Synthesizer 19 zur Verfugung gestellt. Die 
Auswahl entsprechender Steuerprogramme zur Erzeugung eines 
35 Kernspinbildes sowie die Darstellung des erzeugten Kernspin- 
bildes erfolgt liber ein Terminal 21, das eine Tastatur sowie 
einen oder mehrere Bildschirme umfafit. 



30 
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Die vorliegende Erfindung besteht in der Entwicklung bzw. im 
Design neuer Steady-State-Sequenzen durch die eine diffusi- 
onsgewichtete MRT-Bildgebung hinsichtlich der kompli zierten 
5 T 2 -Einwirkung, sowie der Artefakte die durch Bewegung entste- 
hen, wesentlich verbessert wird. Die erfindungsgemaUen Se- 
quenzen werden im Synthesizer 19 gespeichert und von der Se- 
quenzsteuerung 18 durchgef uhrt . 

10 Zum besseren Verstandnis wird im Folgenden anhand von schema- 
tischen Zeichnungen, insbesondere im Zusammenhang mit einem 
Spin-Echo-Experiment, die Theorie und Vorgehensweise der dif- 
fusionsgewichteten Bildgebung erlautert. 

15 Die Aufnahme der Daten in der MRT erfolgt im sogenannten k- 
Raum (Synonym: Frequenzraum) . Das MRT-Bild im sogenannten 
Bildraum ist mittels Fourier-Transformation mit den MRT-Daten 
im k-Raum verkniipft. Die Ortskodierung des Objektes, welche 
den k-Raum aufspannt, erfolgt mittels Gradienten in alien 

20 drei Raumrichtungen . Man unterscheidet dabei die Schichtse- 
lektion (legt eine Auf nahmeschicht im Objekt fest, liblicher- 
weise die z-Achse) , die Frequenzkodierung (legt eine Richtung 
in der Schicht fest, ublicherweise die x-Achse) und die Pha- 
senkodierung (bestimmt die zweite Dimension innerhalb der 

25 Schicht, ublicherweise die y-Achse) . 

* 

V Es wird also zunachst selektiv eine Schicht beispielsweise in 

z-Richtung angeregt. Die Kodierung der Ortsinf ormation in der 
Schicht erfolgt durch eine kombinierte Phasen- und Frequenz- 
30 kodierung mittels dieser beiden bereits erwahnten orthogona- 
len Gradientenf elder die bei dem Beispiel einer in z-Richtung 
angeregten Schicht durch die ebenfalls bereits genannten Gra- 
dientenspulen in x- und y-Richtung erzeugt werden. 



35 



Ein Beispiel fur die Aufnahme von Daten in einem MRT-Experi- 
ment ist in den Figuren 2a und 2b dargestellt. Die verwendete 
Sequenz ist eine Spin-Echo-Sequenz . Bei dieser wird durch ei~ 
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12 

nen 90°Anregungsimpuls die Magnetisierung der Spins in die x- 
y-Ebene geklappt. Im Laufe der Zeit (1/2 T E ; T E ist die Echo- 
zeit) kommt es zu einer Dephasierung A<|) der Magnetisierungs- 
anteile, die gemeinsam die Quermagnetisierung in der x-y- 
Ebene M xy bilden. Nach einer gewissen Zeit (z.B. 1/2 T E ) wird 
ein 180°-Impuls in der x-y-Ebene so einges trahlt , dass die 
dephasierten Magnetisierungskomponenten gespiegelt werden oh- 
ne dass Prazessionsrichtung und Prazessionsgeschwindigkeit 
der einzelnen Magnetisierungsanteile verandert werden, Nach 
einer weiteren Zeitdauer 1/2 T E zeigen die Magnetisierungs- 
komponenten wieder in die gleiche Richtung, d.h. es kommt zu 
einer als „Rephasierung* bezeichneten Regeneration der Quer- 
^jjj^ magnetisierung. Die vollstandige Regeneration der Quermagne- 
tisierung wird als Spin-Echo bezeichnet. 

15 

Urn eine ganze Schicht des zu untersuchenden Objektes zu mes- 
sen, wird die Bildgebungssequenz N-mal fur verschiedene Werte 
des Phasenkodiergradienten z.B. G Y wiederholt, wobei die Fre- 
quenz des Kernresonanzsignals (Spin-Echo-Signals) bei jedem 

20 Sequenzdurchgang durch den At-getakteten Analog-Digital- 

Wandler ADC (engl.: Analog Digital Converter) N-mal in aqui- 
distanten Zeitschri tten At in Anwesenheit des Auslesegra- 
dienten G x abgetastet, digitalisiert und abgespeichert wird. 
Auf diese Weise erhalt man gemali Figur 2b eine Zeile fur Zei- 

25 le erstellte Zahlenmatrix (Matrix im k-Raum bzw. k-Matrix) 
& mit NxN Datenpunkten. Eine symmetrische Matrix mit NxN Punk- 
l ten ist nur ein Beispiel, es konnen auch asymmetrische Matri- 

zen bzw. andere k-Raum-Belegungen erzeugt werden. Aus derar- 
tigen Datensatzen im k-Raum konnen durch Fouriertransf ormati- 

30 on unmittelbar MR-Bilder der betrachteten Schicht mit einer 
Auflosung von NxN Pixeln rekonst ruiert werden. 

Das Auslesen muii in einer Zeit abgeschlossen sein, die klein 
ist gegenuber dem Zerfall der Quermagnetisierung T 2 . Ansons- 
35 ten waren namlich die verschiedenen Zeilen der k-Matrix ent- 
sprechend der Reihenfolge ihrer Erfassung unterschiedlich ge- 
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wichtet: bestimmte Ortsf requenzen wurden liber-, andere dage- 
gen wurden unterbetont werden. 

In einem moglichen Beispiel eines dif f usionsgewichteten Spin- 
Echo-Experimentes wird gemaii Figur 3 ein konstanter starker 
Dif fusionsgradient 25 angelegt. Wird der Spin durch den 90°- 
Puls in die x-y-Ebene gekippt, so kommt er anfangs auf genau 
einer Achse (z.B. der x-Achse) zu liegen. Der Spin beginnt 
nun in dem starken Zusatzfeld des Dif f usions-Gradienten 25 zu 
Dephasieren, d.h. nach einer gewissen Zeit hat er einen ge- 
wissen Winkel § eingenornmen. Dieser zeitliche Winkelverlauf (() 
- auch als Phase (j) bezeichnet - ist gemaii Figur 3 abhangig 
^ von der Dauer (im Falle eines zeitlich begrenzten Diffusions- 
' Gradient-Pulses) und der Starke des anliegenden Diffusions- 

15 Gradienten 25. Nach einer gewissen Zeit (TE/2) wird der 180°- 
Puls eingestrahlt . Zu diesem Zeitpunkt hat die Magnetisierung 
eine Dephasierung von (j)i erreicht. Durch den 180°-Puls wird 
die Magnetisierung in der x-y-Ebene invertiert (-(J>i) und be- 
ginnt zu Rephasieren bis nach der Zeit TE ein maximal mogli- 
20 ches Echo 23 zustande kommt. 

Anders verhalt es sich bei einer dif f usionsbedingten Ortsan- 
derung Kernspin-resonanter Materie (z.B. Wasser-Molekule) 
zwischen Einstrahlung des 90°-Pulses und des 180°-Pulses zu 

25 einem beliebigen Zeitpunkt t. Bewegen sich die Molekule bei- 
spielsweise an einen Ort hoherer Dif f usionsgradientenstarke, 
( so erfahren sie in dem hoheren lokalen Feld eine starkere 

Dephasierung im Vergleich zu der nichtbewegten Materie (Phase. 
4>i) . Verbleibt im vereinf achten Fall die bewegte diffundierte 

30 Materie an dem neuen Ort, so erfahrt deren Magnetisierung 

durch den 180°-Puls eine Invertierung bei einer grofieren Pha- 
se (|)2. Aufgrund des hoheren lokalen Feldes hat die Phase auch 
nach der Invertierung einen zeitlich gesehen steileren Ver- 
lauf, was dazu flihrt, dass die Magnetisierung der diffundier- 

35 ten Materie zu dem Auslesezeitpunkt TE eine Zusatzphase Acf) 
aufweist. Diese hat zur Folge, dass die Quermagnetisierung 
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der bewegten Materie im Ganzen ein geringeres Echosignal 24 
liefert als die nicht-bewegte Marterie. 

Dies kann anhand der Figur 4 weiter veranschaulicht werden: 
5 Auf der linken Seite ist die sich aus Magnet isierungsanteilen 
zusammensetzende Quermagnetisierung der nicht-bewegten Kern- 
resonanten Marterie zum Echo-Zeitpunkt (TE) vektoriell darge- 
stellt. Wie man sieht erreicht die Magnetisierung wieder ih- 
ren Maximalwert, der dem Ausgangswert zum Zeitpunkt des 90°- 

10 Pulses entspricht (Zt=I max ) . Auf der rechten Seite sieht man, 
dass die Magnetisierungsanteile der bewegten (dif f undierten) 
kernresonanten Materie zum Echozeitpunkt unterschiedliche 

^ Phasen A(J> besitzen und die Vektorsumme daher eine geringere 
Magnetisierung ergibt als sie theoretisch moglich ware Zt<I- 

1 5 m ax ) • 



Zusammenf assend kann gesagt werden: Je grower die Diffusion, 
d.h. die Ortsveranderung der kernresonanten Materie im Gra- 
dientenfeld, desto geringer ist das Kernresonanzsignal in 
20 diesem Bereich. Bereiche starker Diffusion machen sich daher 
in einer Signalabschwachung bemerkbar. 

Die dif fusionsgewichtete Bildgebung in der MRT erfordert aus- 
gesprochen hohe Dif f usionsgradientenf eldstarken, da der oben 
25 beschriebene Effekt sehr klein ist. Bei einem zwar starken 

aber konstanten Dif f usions-Gradienten - wie er in Figur 3 ab- 
I gebildet ist - hat man den Nachteil, dass auch das Auslesen 

des Echos unter einem starken Gradienten erfolgt und somit 
gemaJi der Beziehung 



30 



Bandbreite^ 1 



Auslesezeit 



35 



durch eine kurze Auslesezeit - bedingt durch den starken Gra- 
dienten - das Signal mit hoher Bandbreite gemessen wird. Eine 
hohe Bandbreite bedeutet, dass das eigentliche Kernresonanz- 
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signal auf einem sehr breitem Frequenzband gemessen und damit 
sehr viel Rauschen mitgemessen wird. 

Aus diesem Grund sind Dif f usionskodierung und Auslesen des 
5 kodierten Kernresonanzsignals zu trennen. 

Eine dif f usionsgewichtete Spin-Echo-Sequenz-Messung rait ge- 
pulsten Dif fusionsgradienten wurde erstmals von Stejskal und 
Tanner verof f entlicht (Stejskal and Tanner, Journal of Chemi- 
10 cal Physics, 42, 288 (1965)) und ist in Figur 5 schematised 
dargestellt : 

Die pulsformigen Dif fusionsgradienten G D if f haben eine relativ 
zu den anderen Feldgradienten (G S ,G P ,G R ) urn ein vielf aches ho- 
15 here Amplitude (G D iff « lOxGR) , weisen eine definierte zeitli- 
che Breite 5 auf und haben zueinander einen definierten zeit- 
lichen Abstand A der auch als Dif f usionszeit bezeichnet 
wird. 

20 Die dargestellte Sequenz hat zuni einen den Vorteil, dass man 
das Echo mit sehr kleiner Bandbreite auslesen kann, indem ein 
sehr kleiner Auslesegradient verwendet wird. Zum andern legen 
die definierten GroBen G D iff, 5 und A insbesondere die zeitli- 
che Charakteristik der Dif fusions-Messung fest, so dass der 

25 Einfluft der Gradientenpulse bzw. deren Abfolge in einem Dif- 
fusions-gewichteten MRT-Experiment mathematisch durch den so- 
{ genannten b-Wert beschrieben werden kann: 



30 



mm 2 



Dabei stellt y das gyromagnetische Verhaltnis der betrachte- 
ten Spinspezies ( im allgemeinen Wasserstoff ) dar. Der die 
Diffusion beschreibende physikalische Wert, der mit b gewich- 
tet wird, wird allgemein als "Dif fusions-Koef f izient D" be- 
35 zeichnet und stellt gemaii der Beziehung 
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D 



tj a 



mm 



eine Naturkonstante dar, wobei T die Temperatur, r\ die Visko- 
sitat und a der Molekulradius der zu messenden dif f undieren- 
5 den Materie ist. Die genaue Einheit ergibt sich unter Beriick- 
sichtigung eines Proportionalitatsf aktors der hier nicht an- 
gegeben ist. In der Realitat ist eine Dif fusionsmessung durch 
weitere Flusseffekte wie Maschineninstabilitaten, Flussbewe- 
gungsartef akte, usw. kontaminiert , die in dem zu messenden 

10 Dif fusionskoef f izienten D berucksichtigt werden mussen. Aus 

diesem Grund hat sich die Bezeichnung "Scheinbarer Diffusions 

^ Koeffizient" (engl.: Apparent Diffusion Coefficient, ADC) D ADC 
durchgesetzt, der - urn ein dif fusionsgewichtetes MRT-Bild zu 
erhalten - in dif f usionsgewichteten MRT-Experimenten zu mes- 

15 sen ist. 

In der sogenannten Zweipunkt-Methode erhalt man ein diffusi- 
onsgewichtetes MRT-Bild einer Schicht dadurch, dass man zu- 
nachst zwei dif f usionsgewichtete Bilder der gleichen Schicht 
20 mit unterschiedlichen b-Werten (bi und b 2 ) misst, aus diesen 
beiden Bildern pixelweise den D AD c-Wert berechnet und diesen 
in einer Bild-Matrix dieser Schicht darstellt. 

. Das Kernresonanzsignal einer dif f usionsgewichteten Messung 
Wb lafit sich mathematisch darstellen als 

S Diff = S 0 * e'^ * . 

Der Faktor e ( ~ bD) wird allgemein als "Dif f usionswichtung" be- 
30 zeichnet und charakter isiert die durch Diffusion bedingte 

(-^) 

Signaldampfung. Durch den zweiten Dampf ungsf aktor e Tl wird 
der durch T2 beschriebene Zerfall der Quermagnetisierung wah- 
rend der Echozeit TE berucksichtigt. S 0 beschreibt MefJ- und 
Gewebeparameter wie beispielsweise Spindichte, Repetitions- 
35 zeit, Flipwinkel, usw . 
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Durch Quotientenbildung lasst sich der D AD c-Wert jedes Pixels 
allein durch die beiden Meflwerte S x und S 2 jedes Pixels sowie 
durch die beiden unterschiedlichen b-Werte bi und b 2 wie 
5 folgt berechnen: 



Fullt man die Schichtmatrix mit den so berechneten D ADC - 

10 Werten, so erhalt man ein dif f usionsgewichtetes Bild dieser 
Schicht. Unterschiedliche b-Werte gewinnt man iiblicherweise 

^ durch Variation der Dif f usionsgradienten-Parameter G, A und 
5. Die Werte S x und S 2 werden aufgrund der Schnelligkeit des 
Messverfahrens vorteilshaf terweise mit Einzel-Schul5-Echo- 

15 Planarer-Bildgebung (engl. : Single-Shot-Echo-Planar-Imaging, 
SSEPI) gewonnen. Dadurch lassen sich bewegungsbedingte Bild- 
artefakte vermeiden. Allerdings besitzt die SSEPI-Sequenz - 
wie eingangs erwahnt - eine durch die Art der Phasenkodierung 
bedingte starke T 2 - bzw. Phasen- Empf indlichkeit die in ex- 

20 trem starken Bildausloschungs- bzw. Ver zerrungsartef akten re- 
sultiert. Dies muii mangels Alternative derzeit in Kauf genom- 
men werden. 

1989 haben M.H.Cho and C . H . Cho in dem Konf erenzband "Society 
of Magnetic Resonance in Medicine p. 911 Amsterdam 1989" theo- 
retisch gezeigt, dass in der Dif fusionsgewichteten MRT- 
Bildgebung bei Verwendung eines bipolaren Dif f usionsgradien- 
tenpuls und der Kenntnis zweier dif fusionsgewichteter SSFP 
(engl.: Steady-State-Free-Precession) Signale S + und S" - je- 
30 weils gemessen mit und ohne Dif f usionsgradient - eine einfa- 
che Formel zur Berechnung des Dif fusionskoef f izienten D ADC an- 
gegeben werden kann: 



D 



ADC 



= ln(-^-) / (b 2 -b x ) 



D 



ADC 



2*b, 



bip 



In 




wobei 



35 
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S~ das Signal einer PSIF-Sequenz ohne Dif f usionsgradienten, 
Sq das Signal einer FISP-Sequenz ohne Dif f usionsgradienten, 
S~ ijr das Signal einer PSIF-Sequenz mit bipolaren Diffusions- 
gradienten G D ±ff, 

5 S^p das Signal einer FISP-Sequenz mit bipolaren Diffusions- 
gradienten G D ±ff und 

b bip der b-Wert der obigen dif f usionsgewichteten bipolaren 
PSIF- bzw. FISP-Messung ist. 

10 Anschaulich wird also der Quotient aus unterschiedlich er- 

zeugten Bildern gebildet, logarithmiert und mit dem entspre- 
chenden b-Wert b b i P gewichtet. 

\ 

Oben genannte Formel zur Bestimmung des D ADC -Koef f izienten hat 
15 den erheblichen Vorteil, dass die Verwendung von Steady- 

State-Signalen (diese werden mittels der herkommlichen Spin- 
Warp-Phasenkodiertechnik erzeugt) bezuglich der Artefakte, 
unter der eine Dif f usionsgewichtete MRT-Messung mit EPI- 
Sequenzen leidet, unempf indlich ist. 

20 

Die Verwendung eines bipolaren Dif f usionsgradienten stellt 
technisch kein Problem dar; der entsprechende b-Wert b bi der 

jeweiligen Geratekonf iguration fur eine bipolare diffusions- 
gewichtete MRT-Messung kann ohne weiteres angegeben werden. 

2 5 

W Die experimentelle Ermittlung der jeweiligen SSFP-Werte 
(S~,SQ,Sp iJT und Sp iJT ) ist theoretisch durch die zeitliche 

Aufeinanderfolge einer FISP-Sequenz und einer PSIF-Sequenz 
bzw. umgekehrt zu Realisieren. Eine solche Vorgehensweise der 

30 Ankopplung beider MeBsequenzen bedeutet im wesentlichen eine 
Verdoppelung der Mefizeit mit dem damit verbundenen schwerwie- 
genden Nachteil, dass aufgrund der aufterordentlichen Empfind- 
lichkeit von FISP- und PS I F-Sequenzen gegenuber Flufi- und 
Pulsstorungen Artefakte mitgemessen werden die das diffusi- 

35 onsgewichtete Bild in seiner Qualitat erheblich beeintrachti- 
gen. 
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In der vorliegenden Erfindung wird vorgeschlagen, zur Messung 
der fur die obige Formel erf orderlichen SSFP-Signale, beide 
Signale S + und S" simultan zu akquirieren. Eine derartige Se- 
quenz wird als DESS-Sequenz (engl.: Double Echo Steady State) 
5 bezeichnet. DESS verknupft die FISP-Sequenz mit dem zusatzli- 
chen Signalgewinn der PSIF-Sequenz und liefert insbesondere 
fur Gewebe mit langer Relaxationszeit T 2 - wie es bei Flus- 
sigkeitsansammlungen (beispielsweise bei Gelenksver letzungen) 
der Fall ist - einen guten T 2 -Kontrast. Dazu werden ublicher- 
10 weise in derselben Sequenz FISP- und PSIF-Signale erzeugt, 
die anschlieflend addiert werden. 



^ Eine DESS-Sequenz ist in Figur 6 dargestellt. Das FISP-Echo 
wird nach jedem Anregungspuls a durch den ersten rephasie- 

15 renden Puis des Auslesegradienten erzeugt. Insofern stellt 
FISP eine Gradientenechosequenz dar. Der dephasierende und 
rephasierende Phasenverlauf des Quermagnetisierungsvektors 
der FISP-Anregung ist im unteren Teil der Figur 6 hell- 
schraf f iert dargestellt . 

20 

PSIF ist eine riickwarts ablaufende FISP-Sequenz (daher das 
Akronym "PSIF") . Bei der PSIF-Technik wird ausgeniitzt, dass 
jeder HF-Puls nicht nur eine anregende, sondern auch eine re- 
fok^ssierende Komponente besitzt. Betrachtet man mehrere An- 
25 regungen einer PSIF-Pulsf olge, so wirkt die Anregung einer 




Messung ref okussierend auf die Transversalkomponente (Quer- 



magnetisierung) der vorhergehenden Messung: Der dephasierende 
und rephasierende Phasenverlauf des Quermagnetisierungsvek- 
tors der PSIF-Anregung ist in der Figur 6 dunkel-schraf f iert 

30 dargestellt. Es wird ein Echo erzeugt, das mit dem zweiten 

rephasierenden Puis des Auslesegradienten ausgelesen wird. Da 
das PSIF-Signal, das mit der ersten Anregung erzeugt wird, 
erst am Ende des Zyklus zur zweiten Anregung ausgelesen wird, 
ist die Echozeit gleich der doppelten Repetitionszeit . Die 

35 resultierenden PSIF-Echos besitzen daher eine ausgesprochen 
hohe T 2 -Wichtung. Streng genommen fallt PSIF nicht in die 
Klasse der Gradientenechosequenzen; es handelt sich hier 
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vielmehr um eine Spinechosequenz . PSIF ist demnach unempfind- 
lich gegenuber existierenden Suszeptibilitaten . 

Die vorliegende Erfindung besteht darin, die DESS-Sequenz in 
5 der dif fusionsgewichteten MRT-Bildgebung einzusetzen. Dazu 
wird die DESS-Sequenz er f indungsgemaft einmal mit bipolarem 
Dif fusionsgradienten ein zweites mal ohne Di f f usionsgradien- 
ten gemessen, wobei erf indungsgemaii der bipolare Diffusions- 
gradient unmittelbar zwischen dem Auslesen des FISP-Signales 
10 und dem Auslesen des PSIF-Signales geschaltet wird. Die in 

der jeweiligen Sequenz gemessenen FISP- und PSIF-Signale (S + 
und S") werden erf indungsgemal} nicht addiert sondern sepa- 
^ riert gespeichert (beispielsweise in dem Anlagenrechner 18) 
V und in die obige Formel eingesetzt. Der ebenfalls in der For- 
15 mel benotigte b-Wert b bip ist durch die Amplitude und Puls- 

breite des bipolaren Dif fusionsgradienten gegeben. Hierdurch 
kann fur jedes Pixel der gemessenen Schicht der D ADC -Wert be- 
rechnet werden. Die Matrix-Darstellung aller D ADC -Werte der 
Schicht ergibt letztendlich ein artef akt f reies dif f usionsge- 
20 wichtetes Bild. 

Eine erste Ausf uhrungsf orm der erf indungsgemafien Sequenz, in 
der ein FISP- und ein PSIF-Signal vor bzw. nach einem bipola- 
ren Dif fusionsgradientenpuls gemessen wird, ist in Figur 7 

25 dargestellt. Wahrend des HF-Anregungspulses a wird der 

p Schichtselektionsgradient GS geschaltet. Unmittelbar daran 

schlieftt sich nach einem dephasierenden Puis des Auslesegra- 
dienten G A - bei gleichzeitiger Phasenkodierung durch den 
Phasenkodiergradienten G p - ein rephasierender Auslesepuls 

30 an, der ein FISP-Echosignal erzeugt. Die Zeit zwischen HF- 

Puls und FISP-Echo wird als Echozeit TE+ bezeichnet. Nach der 
Akquirierung des FISP-Echos wird der Dif fusionsgradienten- 
Puls mit der Gradientenamplitude G 0 und der Pulsbreite 5 bi- 
polar geschaltet und daran anschlieflend durch einen weiteren 

35 Gradientenpuls des Auslesegradienten das PSIF-Echo S* ausge- 
lesen. Noch vor dem nachsten HF-Anregungspuls wird der PSIF- 
Auslesegradient rephasierend geschaltet, wahrend die Phasen- 
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kodierung des Phasenkodiergradienten durch einen entsprechend 
invertierten Phasenkodiergradientenpuls aufgehoben wird. Die 
Zeit zwischen Akquirierung des PSIF-Echos und dem nachsten 
HF-Puls wird als te- bezeichnet, so dass die Echozeit des 
5 PSIF-Echos TE- durch TE-=2*TR-te angegeben werden kann. Dabei 
ist TR wie ublich die Repetitionszeit , d.h. die Zeit zwischen 
auf einanderf olgenden HF-Anregungspulsen . 

Streng genommen setzt sich ein PSIF-Echo aus der Refokussie- 

10 rung mehrerer nicht unbedingt auf einanderf olgender HF-Pulse 
zusammen, was hinsichtlich des Phasenverlauf es sowohl der 
Quer- als auch der Langs-Magnetisierung zu einer komplizier- 

^ ten Auffacherung (engl.: Splitting) f uhrt . Diese Auffacherung 
kann vereinfacht durch sogenannte Echopfade dargestellt wer- 

15 den, deren jeder einen anderen Zeitverlauf, insbesondere eine 
andere Dif f usionszeit Ai aufweist. i indiziert die Anzahl der 
HF-Pulse die bis zur Entstehen des jeweiligen PSIF-Echos ei- 
nes bestimmten Echopfades erfolgt sind. Durch die Bipolaritat 
des Dif fusionsgradientenpulses sowie durch das simultane Aus- 

20 lesen beider Echos FISP und PSIF durch mehr oder weniger in- 
verse Auslesegradientenpulse wird die Komplexitat insbesonde- 
re der PSIF-Signalgenerierung kompensiert . Dieser Sachverhalt 
spiegelt sich wider in der Einfachheit der Formel nach Cho . 

25 Vergleicht man die erf indungsgemaJJe Sequenz von Figur 7 mit 

| der herkommlichen DESS-Sequenz in Figur 6, so fallt auf, dass 
die Totalsyxnmetrie der Pulse der DESS-Sequenz in Figur 7 
nicht gegeben ist. Dies mufi nicht so sein, ist aber ein vor- 
teilhafter Aspekt der vorliegenden Erfindung und hat folgen- 

30 den Grund: Das FISP-Echo S + hat von Natur aus ein hoheres 
Signal-zu-Rausch-Verhaltnis (engl. : Signal-to-Noise-Ratio, 
SNR) als das PSIF-Echo S~, welches aufgrund des relativ ge- 
ringen Ref okussieranteils eines HF-Anregepuls verhaltnismaftig 
klein ist. Das Auslesen des PSIF-Echos S" erfordert deshalb 

35 einen zeitlich gesehen langen Auslesepuls danger ADC- , wobei 
hier mit ADC wieder der auslesende Analog-Digital-Converter 
gemeint ist) . Das Auslesen des FISP-Echos S + kann in einer urn 
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ein Vielfaches kurzeren Zeit erfolgen, wobei der Auslesegra- 
dient dafiir eine entsprechend hohere Amplitude aufweisen muii . 
In anderen Worten bedeutet dies, dass das FISP-Echo ntit hoher 
Bandbreite gemessen werden kann - ohne SNR-Verlust, wodurch 
5 insgesamt, uber den gesamten Zeitraum der Messung gesehen, 
die Messzeit erheblich reduziert wird. 

Werden FISP-Echo und PSIF-Echo gemafi der ersten Ausfuhrungs- 
form der vorliegenden Erfindung mit verschiedenen Bandbreiten 

10 gemessen, so kann dies (beispielsweise aufgrund der chemi- 

schen Verschiebung) dazu fuhren, dass das FISP-Bild gegenuber 
dem erzeugten PSIF-Bild ortlich verschoben akquiriert wird. 

|~ Bei der Quotientenbildung nach der Formel von Cho ftihrt dies 

* zu deutlichen und daher storenden Artefakten. Um dies zu ver- 

15 meiden besteht eine zweite Ausf uhrungsf orm der Erfindung dar- 
in, das PSIF-Echo mit einer sogenannten n-Multi-Echof olge (n- 
Multigradientenechozug) auszulesen, wie es in Figur 8 fur n=3 
dargestellt ist. Eine solche n-Multi-Echof olge besteht aus 
einer Anzahl ( typischerweise aus 2n) alternierender Gradien- 

20 tenpulse der gleichen Pulsbreite wie der des FISP-Echo- 

Auslesepulses . Dies liefert n PSIF-Echo-Signale die mit der 
gleichen Bandbreite wie das FISP-Echo gemessen werden. Durch 
Addition und anschliefiende Mittelung (z.B. Quadratsummenbil- 
dung) der einzeln akquirierten PSIF-Echo-Signale ist dennoch 

25 die gleiche Empf indlichkeit (SNR) gegeben wie bei einem lan- 

| geren Auslesen (ADC- ) gemafl Figur 7. 

Ferner ist es moglich, die eben beschriebene erf indungsgemaUe 
simultane bipolar-dif f usionsgewichtete Akquisitions-Technik 

30 gemafi Figur 9 auf das bekannte MRT-Pro j ektions-Rekonstruk- 
tionsverf ahren anzuwenden. Das Pro j ektions-Rekonstruktions- 
verfahren (Lauterbur, New York 1973) ist eines der altesten 
bildgebenden Verfahren in der Magnetresonanztomographie und 
erfahrt derzeit aus unterschiedlichen Grunden eine Renais- 

35 sance. Es ist ein Verfahren, welches die Berechnung eines di- 
gitalen Bildes aus einer Serie von Projektionen ermoglicht. 
Dazu wird (wie bei alien anderen bildgebenden Verfahren in 



200204698 



23 

der MRT, insbesondere in der 2- bzw. 3-dimensionalen Fourier 
Transf ormations-Tomographie) gemafi Figur 9 durch einen 
Schichtselektionsgradienten GS eine axiale Schicht in dem zu 
untersuchenden Volumen angeregt bzw. selektiert. Unmittelbar 
5 danach wird ein gleichbleibender Gradient G B iid urn kleine Win- 
kelinkremente A(|) gedreht, wobei jedes Mai eine Projektion er 
zeugt wird. Die Rekonstruktion eines aus NxN-Pixeln bestehen 
den Bildes erfordert demnach N unabhangige Pro j ektionen, die 
ihrerseits durch N Punkte definiert sind. Die so in Polarko- 
10 ordinaten gewonnenen Pro j ektionen konnen kartesisch durch 

zwei orthogonale Gradienten G x und G y beschrieben werden, wo- 
bei ohne Beschrankung der Allgemeinheit gilt: 

V G y = G md sin(^) , 

15 G x = G md cos(^) . 

Die Kombination der erf indungsgemaBen Sequenz mit dem Projek 
tions-Rekonstruktions-Verf ahren bewirkt eine deutliche Reduk 
tion von Bewegungsartef akten . 

20 

Zum Schlufi soil durch eine grobe Abschatzung die bipolare 
Dif fusionsgewichtete MRT mit der monopolaren verglichen wer- 
den. Der b-Wert der klassischen Stej skal-Tanner-Dif fusions- 
sequenz mit monopolaren Gradientenpulsen (wie sie in Figur 5 
25 dargestellt ist) wird gemaft der bereits weiter oben angefuhr 
ten mathematischen Beziehung 

30 berechnet. Wenn die Dif f usionszeit A gleich der Pulsdauer 5 
des Dif f usionsgradienten ist (A=5) , gilt 

2 

b mono = ^r 2 G D iff 2 ^ • 
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24 

Wird bei der bipolaren Dif f usions-Gradientenschal tung insge- 
samt die gleiche Pulsdauer verwendet, so gilt ferner 

Dies bedeutet fur den b-Wert der erf indungsgemafien bipolaren 
Sequenz 

Kp = 2 mal ~b mono = ^b mono 
Demnach lafit sich der in der Formel nach Cho benotigte b-Wert 



if9r bbip folgendermafien berechnen: 



Obige Abschatzung zeigt, dass fur einen vergleichbaren Gewe- 
bedif fusionsef f ekt mit monopolarer Technik im bipolaren Ver- 
fahren Dif f usionsgradienten mit einer um ein Vielfaches hohe- 
ren Amplitude angelegt werden miissen. 



Die Bipolaritat des Dif f usions-Gradientenpulses G D iff der er- 
f indungsgemaJien Sequenz hat zudem den Vorteil, dass Nichtli- 
nearitaten bzw. Verzerrungen des Dif fusionsgradienten-Pulses 
^ - bedingt durch Wirbelstrome die in metallischen Leitflachen 
25 (Kalteschild, Shimbleche, etc.) durch Auf- und Abbau des Gra- 
dientenf eldes induziert werden - kompensiert werden. 
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Patent ansprilche 

1. Verfahren zur Ermittlung eines dif f usionsgewichteten Bil- 
des in der dif f usionsgewichteten MRT-Bildgebung, 
5 aufweisend die folgenden Schritte: 

a) Messen und Speichern eines nichtdiff usionsgewichteten und 
eines dif f usionsgewichteten Datensatzes mittels einer 
DESS-Sequenz (Double Echo Steady State Sequenz) , wobei fur 
den nichtdif f usionsgewichteten Datensatz zwei Auslesegra- 

10 dienten nacheinander geschaltet werden und fur den diffu- 

sionsgewichteten Datensatz zwischen zwei Auslesegradienten 
ein bipolarer Dif fusionsgradientenpulszug geschaltet wird, 

b) Berechnen eines dif f usionsgewichteten MRT-Bildes auf Basis 
des nichtdif f usionsgewichteten und des dif f usionsgewichte- 

15 ten Datensatzes sowie auf Basis eines die dif f usionsge- 

wichtete Messung charakterisierenden Wertes. 



20 



2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 

dass der bipolare Dif fusionsgradientenpulszug aus einem posi- 
tiven und einem negativen Dif f usionsgradientenpuls mit der 
jeweiligen Amplitude G 0 bzw. -G 0 sowie der gleichen Pulsbrei- 
te 5 besteht, und wobei beide Dif f usionsgradientenpulse un- 
mittelbar hintereinander folgen. 

25 

*4) 3 - Verfahren nach Anspruch 1 oder 2, 

dadurch gekennzeichnet, 
dass die Berechnung des dif f usionsgewichteten MRT-Bildes 
durch Bilden des Quotienten einer Kombination des diffusions- 
30 gewichteten und des nichtdif f usionsgewichteten Datensatzes, 
anschliefiender Logari thmierung des Quotienten sowie Gewich- 
tung mit dem der dif f usionsgewichteten Messung entsprechenden 
b-Wertes erfolgt. 

35 4. Verfahren nach Anspruch 3, 

dadurch gekennzeichnet, 
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dass das dif f usionsgewichtete MRT-Bild durch pixelweise Dar- 
stellung eines ADC-Koef f izienten D ADC gebildet wird der aus 
den akquirierten Datensatzen gemafi der Formel 

pixelweise ermittelt wird, wobei 

Sq und S~ den Datensatz der nichtdif f usionsgewichteten Mes- 
sung in Form von FISP-Echo-Signalen und von PSIF-Echo-Signa- 
10 len darstellen sowie S^ iJf und S& ff den Datensatz der diffusi- 

onsgewichteten Messung in Form von FISP-Echo-Signalen und von 
dp PSIF-Echo-Signalen darstellen, und wobei 

bbip den die dif f usionsgewichtete Messung charakterisierenden 
Wert darstellt der nach der Formel 

15 

durch die Beschaf f enheit des bipolaren Dif f usionsgradienten 
und dem gyromagnetischen Verhaltnis y der betrachteten kern- 
20 resonanten Spinspezies gegeben ist. 

5. Verfahren nach Anspruch 4, 

dadurch gekennzeichnet, 
^ dass das S + -Signal im Vergleich zum S"-Signal mit hoherer 
25 Bandbreite gemessen wird. 

6. Verfahren nach Anspruch 4, 

dadurch gekennzeichnet, 

dass beide Signale S + und S" mit gleicher grolier Bandbreite 
30 gemessen werden, wobei das S~-Signal durch einen Multigra- 

dientenechozug mehrmals gemessen wird und liber die gemessenen 
S"-Signale gemittelt wird. 

7. Verfahren nach Anspruch 6, 

35 dadurch gekennzeichnet, 
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dass die Mittelung durch ein Quadratsummenver f ahren erfolgt. 

8. Verf ahren nach Anspruch 1 bis 7, 
dadurch gekennzeichnet, 

dass die Akquirierung der Datensatze S~ iff , S^ iff , S~ , S* durch 

ein Proj ekt ions -Rekonstrukt ions -Verf ahren erfolgt • 

9. Kernspintomographiegerat zur Durchfuhrung eines Verfahrens 
nach den Anspruchen 1 bis 7. 
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Zusammenf assung 

Verfahren zur Ermittlung des ADC-Koef f izienten in der Diffu- 
sions-gewichteten Magnetresonanz-Bildgebung bei Verwendung 
5 von Steady-State-Sequenzen 

Die vorliegende Erfindung beansprucht ein Verfahren zur Er- 
mittlung eines dif fusionsgewichteten Bildes in der diffusi- 
onsgewichteten MRT-Bildgebung, aufweisend die folgenden 
0 Schritte: 

a) Messen und Speichern eines nichtdif fusionsgewichteten und 
eines dif fusionsgewichteten Datensatzes mittels einer 
DESS-Sequenz (Double Echo Steady State Sequenz) , wobei fur 
den nichtdif fusionsgewichteten Datensatz zwei Auslesegra- 

5 dienten nacheinander geschaltet werden und fur den diffu- 

sionsgewichteten Datensatz zwischen zwei Auslesegradienten 
ein bipolarer Dif f usionsgradientenpuls geschaltet wird, 

b) Berechnen eines dif fusionsgewichteten MRT-Bildes auf Basis 
des nichtdif fusionsgewichteten und des dif f usionsgewichte- 

0 ten Datensatzes sowie auf Basis eines die dif f usionsge- 

wichtete Messung charakterisierenden Wertes. 
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FIG 2A 

90°-lmpuls 180°-lmpuls Echo 




FIG 2B 
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FIG 5 
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FIG 7 
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FIG 8 
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FIG 9 
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